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Abstract 



We present a model that explains two phenomena, recently observed in high- 
mobility Si-MOS structures: (1) the strong enhancement of metallic conduc- 
tion at low temperatures, T < 2 K, and (2) the occurrence of the metal- 
insulator transition in 2D electron system. Both effects are prescribed to the 
spin-orbit interaction anomalously enhanced by the broken inversion symme- 
try of the confining potential well. 
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Introduction. Recently, in experiments with high mobility Si-MOS structures, a strong 
drop in resistivity p(T) has been found [0 as temperature decreases below ~ 2 K. This 
effect is evidently in disagreement with the conventional interpretation of the one-parameter 
scaling theory (OPST) ||, according to which all states in 2D system at zero magnetic 
field should be localized in the limit of T — > 0. The subsequent scaling analysis of the 
temperature and electric field dependencies of the conductivity has revealed a critical 
behavior, typical for the metal-insulator transition. Finally, convincing evidence for the 
existence of the extended states in 2D system at zero field has been obtained in experiments 
in magnetic field, in studies of the quantum Hall effect to insulator transitions |EJ. The 
extended states which in high magnetic field are located in the centre of the corresponding 
Landau bands, at decreasing field were found to remain in a finite energy range, giving rise 
to a mobility edge. 

The experimental results thus suggest the existence of a true metallic state and of the 
metal-insulator (M-I) transition in 2 dimensions. These results are in apparent contradic- 
tion with the conventional OPST, and the origin of the metallic state remains puzzling. 
In this work, both experimental findings are explained as a consequence of the spin-orbit 
interaction enhanced by the broken inversion symmetry. The suggested model provides a 
good agreement with the experimental data on the temperature dependence of the resistivity 

Analysis of the experimental results. Fig. 1 shows a set of the curves p vs T 
typical for the high-mobility samples 0, at different electron densities n s . At T > 2 K, the 
resistivity, p, increases slowly as temperature decreases, the latter is characteristic for the 
weakly localized regime. At lower temperature, p drops sharply for all curves belonging to 
the "metallic" range of densities, n s > n c . 

The resistance drop is observed at densities in the range from n c to ~ 1.4 x n c . The 
critical density, n c is sample dependent and is equal to 9.2 • 10 10 cm -2 for the sample shown 
in Fig. 1. The drop in p(T) diminishes with decreasing sample mobility, and is almost 
replaced by a conventional rise in p at T — > in the sample with 8 times lower mobility, 
p = 5000 cm 2 /Vs. The latter behavior is consistent with that reported in earlier studies on 
low-mobility samples ||. 

Empirical fit of the data. The p(T)- curves in Fig. 1 may be fitted well by an empirical 
dependence which summarizes the scattering probabilities of two processes: 

p(T)=p + Pl exp(-T*/T). (1) 

The first term is independent of temperature, while the second one describes a scattering 
through an energy gap, A = ksT*. The curves shown in Fig. 1 by thick continuous lines 
were obtained using two fitting parameters for each density, pi/po an d T*. 

Possible microscopic mechanisms. As seen from Fig. 1, the characteristic tempera- 
ture, T*, is of the order of 2 K in high mobility samples. Searching for a proper microscopic 
mechanism, we find two small energy gaps intrinsic to Si-MOS structures at H = 0: the 
valley splitting A„ « 2.4 K and the zero-field spin-gap A S (H = 0) ~ 3.6 K ||. 

It seems attractive to link the resistance drop to the transitions between the two electron 
valleys, located close to the X-points in the Brillouin zone 0. However, the intervalley "um- 
klapp" scattering would hardly occur, since it requires a combination of the reciprocal lattice 
vectors of the very high order. On the other hand, the phonon-induced intervalley transition 
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would require participation of high energy phonons, E p h ~ 10 4 K, and is therefore unlikely 
at low temperatures. Electron tunneling, as the intervalley transition mechanism, would not 
lead to a strong temperature dependence of scattering. 

Spin-orbit splitting and interaction effects. In the one-electron approximation, the 
spin-orbit interaction is described by the Hamiltonian J7|: 

^ = 4^ [W(r)XpK (2) 

where p and a are the momentum and spin operators, correspondingly. For the 2D 
electron system in Si, the contribution of the bulk crystal potential in VV^ is small (g* ~ 2), 
and the lack of inversion symmetry of the triangular confining potential U(z) plays the major 
role. 

This lifts the spin- degeneracy at zero magnetic field, and leads to the appearance of a 
linear term in the energy spectrum of 2D electrons ||: 

E±(k) = ^±ak, (3) 

The corresponding spin-gap 

A s = E+ - E~ = 2ak F , (4) 

can be viewed as the difference in energy for the electron states with spin directed in the 
plane but to the left and right side with respect to kp, or, equivalently, along and opposite to 
the effective magnetic field H* ~ (l/m*c) [kpxVC/] in the frame related to electrons moving 
in 2D plane with Fermi velocity. 

We suggest that the empirically determined energy gap, A originates from A s , and is 
equal to pA s (where p ~ 1), whereas the temperature independent contribution to the 
resistivity, po, is related to the spin- independent scattering. Finite quantum relaxation time 
r q and the corresponding level broadening Ti/r q should reduce the effective gap: 

k B T* =pA s -h/r q . (5) 

Comparison of the model with the experimental results. The spin splitting 
A S (H = 0) = 3.6 K for the same Si-MOS structures was determined from the magnitude of 
the quantum oscillations of the chemical potential ||, extrapolated to zero field. It is clear 
from Eqs. (3) and (4), that this value corresponds to E~ in = —2m*a/h 2 . With this result 
we obtain a = 1.79 • Kcm and the total energy spectrum becomes known. 

In order to estimate the quantum level broadening, T = h/r q , the temperature depen- 
dence of the diagonal resistance was measured in the quantum Hall effect regime with Fermi 
energy adjusted to the Zeeman energy gap at v = 6 in the field of H = 0.75 T. As a result, 
we have obtained an estimate, F = 2.8 K for n s = 10.84 ■ 10 10 cm -2 (here we presumed 
the Gaussian level broadening and T to be the full width). Comparing the model effective 
energy gap, A = p A s — T with the empirical value, T* = 0.96 K obtained in the fit at 
n s = 10.84 x 10 10 (the 5th curve from the bottom in Fig. 1), we eventually find p = 0.46 ||. 

The empirical energy gap, k B T*, decreases to zero at n s = n c , as seen in Fig. 1. In the 
above model, Eq. (5), this occurs because (i) the level broadening increases, T oc l/n s and 
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(ii) the spin splitting, A s = 2akp, diminishes oc r^J 2 . The empirical fitting parameter F 
shown in Fig. 2 vs electron density, rises indeed as density decreases. The total fit in Fig. 
1 is in surprisingly good agreement with the experimental data, despite the very simplified 
character of the above model. 

As n s decreases and approaches n c , the resistance drop starts at lower temperatures. 
The weak decrease in p(T) noticeable at T > 4 K is presumably due to the weak localization 
corrections 5p oc — log(T/T ) This effect was ignored in the above model and the p{T) 
points corresponding to the negative dp/dT were not fitted; these points are connected by 
dashed lines in Fig. 1. 

Fig. 2 shows also two relaxation times: r p calculated from the mobility in the T = 
limit, and r q = h/T. It is noteworthy that both r p and r q diminish almost linearly as density 
decreases (but at n s > n c ) and independently of each other. The momentum relaxation time 
t p provides the necessary resistivity value, p ~ h/e 2 at the critical density whereas 
r q provides the effective spin-gap equal to at n s = n c . As a result, t p decays faster and 
becomes smaller than r q with decreasing n s . At the critical density, p is of the order of 
0.1 m 2 /Vs and r ~ 5K, that corresponds to r p = 0.17 x 10~ 12 s, and r q ~ 1.5 x 10~ 12 s. 
The conclusion on interception of the two curves, T q {n s ) and r p (n s ) is model independent, 
whereas the numerical values of r q {n s ) and the interception point depend on the model 
chosen for level broadening. For instance, the interception occurs at 11.5 or 13 ■ 10 10 cm -2 
for Gaussian or Lorenzian broadening correspondingly. 

Metal-Insulator transition at H = 0: Spin-orbit interaction and symmetry 
effects. It is not only the low-temperature resistance drop but the total scaling behavior 
strongly depends on the symmetry. The corresponding universality classes of the symmetry 
for random systems were established by Dyson |TT|. In the presence of the spin-orbit (SO) 
interaction, the orthogonal symmetry of the system is replaced by the symplectic symmetry. 
Correspondingly, the level-repulsion exponent in the random matrix statistics [O changes 
from u = 1 to oj = 4. It appears therefore, that states are less easily localized in systems 
with large to. 

The effect of the spin-orbit interaction on weak localization was studied both theoretically 
and experimentally ||12|| . For the strong localization regime, there have been suggestions [13| 



that M-I transition can occur in the presence of a strong SO interaction. The scaling function 
in 2D was found to behaves asymptotically like (3(G) ~ —a/G in the high conductance limit 
G>1, with a > in the orthogonal and a < in the symplectic case |TJ]]. The j3 -function 



in the symplectic case may thus become positive at sufficiently large G. As disorder increases 
and conduction G decreases, all states will be localized even in the symplectic case. The 
critical level of disorder and the critical conduction G c correspond to the point at which 
P(G C ) =0. 

The behavior of the symplectic /3(G)— function in 2D is qualitatively consistent with the 
experimental data presented in Fig. 1, in the vicinity of n c . As temperature, or broadening 
increases, the energy relaxation time r e appears as a cut-off parameter and l/r e may become 
larger than the inverse spin relaxation time l/r e ^> l/r s . Then the system would again 
behave as in the orthogonal symmetry case. The natural measure of the SO interaction 
strength is the spin-orbit gap A s given by Eqs. (4) and (5), whereas as an estimate for 
disorder we adopted h/r s ~ ^/r g = V. Thus, one may expect the M-I transition would 
manifest in those samples where A s /r > 1, which is also consistent with occurrence of the 
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transition in the samples with peak mobility larger than 5000 cm 2 /Vs 

Discussion. The above model explains why the resistance drop is seen only in low- 
disordered samples with large r q - values ||15|| . The effect is dependent also on the symmetry 
of the potential well. This provides a key for testing the driving mechanism. As for other 
systems, the zero field spin-gap in GaAs/Al(Ga)As is smaller by a factor of 10 -r 100, due 
to the smaller g*-factor value and much smaller VU [|16|]. Thus, even in ideal samples with 
zero broadening, the resistance drop may occur at temperatures 10 100 x lower than those 
for the Si-MOS structures. In accord with this, no signatures of the resistance drop were 
revealed in recent measurements on GaAs/Al(Ga)As heteroj unctions at temperatures down 
to 20mK |7J. 

Recently, there have been suggestions on other possible collective mechanisms, such 
as Coulomb interaction flEfl , spin-triplet pairing and non- Fermi-liquid behavior [50 



However, the corresponding models are not developed yet to provide a comparison with the 
experimental data. 

Summary. It seems likely that the recently observed metal-insulator transition in high- 
mobility Si-MOS structures is the first experimental manifestation of the spin-orbit interac- 
tion induced transition in 2D. The enhancement of the metallic conduction in these samples 
at low temperatures fits the same framework. Strong SO interaction energy relative to the 
level broadening, and broken inversion symmetry are favorable for the 2D metallic state. 
The Coulomb interaction in this model provides the small level broadening at density down 
to n c . In the recent experiments [^l],[2l|, the 2D metallic phase was found to be easily de- 
stroyed by the in-plane magnetic field; this is a strong evidence for the spin-related origin 
of the 2D metal. 
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FIGURES 



FIG. 1. Typical temperature dependencies of the resistivity for a high- mobility sample JT],||]. 
Electron density, from bottom to the top, is equal to 13.69, 12.81, 12.15, 11.50, 10.84, 10.18, 9.53 
8.87, 8.21, 7.55, 7.12 xlO 10 cm" 2 Q[§. Solid lines " simulations, as discussed in the text. 

FIG. 2. Level broadening, T, and relaxation times, as determined from the fit, Fig. 1, in the 
vicinity of n c . 
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